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ABSTRACT 

• • . • I . 

A design for a Mars Roving Vehicle has been underway at 
Rensselaer Polytechnic Institute for the past four years. The basic 
configuration has been developed. It now remains to design and 
develop the deployment mechanisms, and then to optimize the total 
design. 

A basic mathematical model of the vehicle Is a valuable tool 
in the optimization of the design. This project report presents such 
a three dimensional model considering three degrees of freedom. 

In addition, the physical characteristics of the 0.4 scale RPI-MRV 
are presented along with the basic dynamic responses. 
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CHAPTER 1 * 

INTRODUCTION 

Although the tnajor point of emphasis of today's outer space 
program has been shifted during the past year or two, the intriguing 
question of the possibility of some form of life on the planet Mars 
remains. A NASA unmanned mission to our Earth's nearest neighbor is 
currently planned for approximately the 1980 time frame. The basic 
objective of the mission is to deploy scientific packages and to 
conduct scientific experiments. 

It was learned from the Viking Mission that the useful capability 
of a stationary lander is severely limited. Because of topographic 
requirements, a landing site is not necessarily the most interesting 
from a scientific viewpoint. In addition, because of the terrain 
disturbance caused by the landing maneuver, it is desirable to conduct 
scientific explorations in undisturbed areas remote from the landing 
site. Therefore, it is desirable to place a mobile package on the 
surface of Mars. For this purpose, RPI's Martian Roving Vehicle (MRV) 
is being developed (Plate 1). ^ 

A current NASA sponsored project has as its objective the design 
and development of the complete MRV including the locomotion and 
navigation systems. The subject of this paper is the dynamic analysis 
of the RPI design of the MRV. 

The dynamic analysis of the MRV was approached by two methods: 

(1) mathematical modeling 

(2) physical testing 

Initially the equations of motion of a three dimensional, three- 
degree -of- freedom (rol 1 , pitch, vertical displacement) model are 
established. Then, the measurement of the physical characteristics 
of the 0.4 MRV model and dynamic tests were conducted. Finally, a 




RPI-MRV 0.4 Scale Model 




comparison of actual response with the mathematically predicted 
response is made. The result shows the mathematical model to be a 
basic design tool from which gross vehicle motions may be predicted. 
It also may be expanded upon to develop a more sophisticated mathe- 
matical model. 
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• ■ CHAPTER 2 ■ 

THE MATHEMATICAL MODEL 

The 0,4 scale model of the RPI-MRV is shown in Plate 1. Elasticity 
and damping in the suspension of the vehicle is provided almost entirely 
by the, all metal, toroidal, elastic wheel. The rear wheel struts are 
connected to a torsion bar, and the front wheels are attached to a 
flexible front axle. Since the stiffness of the rear torsion bar 
and the front axle is much greater in comparison to the wheel stiff- 
ness, only the wheel elasticity is considered in the mathematical 
model presented herein. ' Figure 1 shows the schematic model. The wheels 
have been replaced by springs and dashpots. This is the schematic from 
which the mathematical model is developed. 

In the development of an initial approach, a highly useable 
mathematical model is the primary objective. Since coupling between 
modes of motion does exist and is a problem in the physical model, 

the idea of a planar mathematical model was discarded in favor of. 

( 

one considering three dimensions, • 

Again for simplicity, three degrees of freedom are considered. 
Vertical displacement is an obvious choice since the MRV bouncing 
across Martian terrain will be carrying a scientific payload. Be- 
cause of the "dragster" design of the RPI-MRV, the vehicle center 
of gravity is located near the rear wheel base far from the center of 
suspension, A natural pitching motion is therefore introduced by 
the vehicle design. As, a result, the mathematical model also 
includes pitch motion. Finally, the movement of the MRV over rocks 




and potholes easily introduces a rolling motion. Therefore, with the 
three degrees of freedom chosen to be vertical displacement, roll, 
and pitch rotation, the development of the equations of motion 
. proceeded. 

The Inertia Forces . 

In order to greatly simplify the equations of motion, the MRV 
is considered to be a two-mass system: a rolling (sprung) mass, 

and a non-rolling (unsprung) or fixed mass. Additionally, if a 
moving coordinate system, x'y'z' fixed to the vehicle is adopted, the 
inertial terms in the equations of motion become much easier to 
determine. Therefore, the coordinate system of Figure 2 is used. 

The xyz axis system is fixed to the unsprung portion of the vehicle. 

The vertical y axis is located by the c.g. of the complete 
vehicle which, for the MRV, is approximated by the location of the 
rolling mass (since m^>m^). The x axis is the roll axis of the 
vehicle, located by the kinematic properties of the suspension. The 
unsprung mass is assumed to be located on this axis. The origin of 
the x'y'z' axes coincides with the sprung (rolling) mass. 

The inertia forces written about the x'y'z' axes are as follows: 

^ Fy* - 

Assuming: (1) Near symmetry in the x'z' plane, I , ,=0 

X z 

(2) Ix'y' 

(3) Ix'y*“^xy^^ since the y axis is chosen through the c.g 
equations (1) become: 
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(to,,) (2) 

For small displacements, the roll and pitch velocities about the 

inclined x'y'z' axes will be respectively equal to the roll and pitch 

* 

velocities about the xyz. axes (i.e. . This approxima- 
tion along with the following expressions: 

• ’ ■ . 


enable equations (2) to be translated to the xyz axes: 


SFy - mf.y 

^ Mjj = 


(3) 


With a change of notation using 
^ZT^Z, 


(4) 


The final inertia forces take the form: 

^'^x“^*x®x (5) 


The External Forces » 

The weight of the sprung mass merely causes a steady static 
deflection of the wheels about which the dynamic oscillations take 


place. Therefore the external forces of the idealized schematic model 



result wholly from the spring and damping forces of the wheel. 

These forces are considered linear. As indicated later, this assump- 
tion proves to be valid. 

To determine the external forces and moments, the vehicle model 
assumes two sets of virtual displacements: 

(1) roll with vertical displacement 

(2) pitch with vertical displacement 

Assuming a roll with vertical displacement (Figure 3), the vertical 
forces become: 

d( ) 

where s denotes the operator . The displacement of the front 
wheel axle is independent of the roll angle displacement because of 
the pin- joint connection between the frame and the front axle. This 
accounts for the last term in equation (6) not being dependent uponG . 

A 

The moments resulting from the displacement of Figure (3) are: 

“(J?j<-h©x)O;,<'C,sXx+i’j0x) + vv»-{^>--/)©x .(7) 

Again, as a result of the pin- joint mentioned above, the front 
suspension imparts no rolling motion to the vehicle. Using the small 
displacement approximations, and, h>>y, equations 

(6) and (7) become: 

Fy = -(lf,+C,s)[(y-^,©J+{y+i^©J+3y] (8) 

Mx= (k,TC,s)[y(J?,-J?,)-ex(J?®-»Jp/}J+Wrk©^ (9) 

Figure 4 shows a pitch angle with vertical displacement. The 




resulting force and moment equations become: 

Fy= "J?(k,+C,s)[(y-a©J + (y-+b6^)J ppj 

Mz~ *Ct|-*C|s) Ua-h©£)(y-c«©J-Cl>+b0e)(y+be!2^ 

*-Wr(h-y)&x. (11) 


Again with the small displacement approximations, b>?-h9^, a»h0^ 
equations (10) and (11) simplify to: 

Fy ^ ^ (|r, + +(1, 

The Complete Equations . 

The complete equations of motion for the three-dimensional, 
three-degree -of- freedom model result from the combination of 
equations (5), (8), (9), (12), and (13): 

M/Y ~ "^*‘i+C,s)[(y-J,G^) + {ytjPj©,)-tJ©^ (t-a)+'/y] 

^xx 4 = { k 1+ c,s) [y " ©X -f/)] + iVr Ji ©X 
^(fet+CiS}^ y (a-b) - Si(a^+fc^J + VV(-h©i 

These equations may be simplified as follows: 

y -t Q„y + Q,J, y + b„ ©X + t>„©x + C„ ©a + c,3©i - O 

^ + b„ex+ Q« y + q,, y = O 

©i+ c^, 6^. + (^©2 +Qj, y 4 y = O 

where 


ail=-4ci/nir 

bii=(52-J!pc^/m^ 

ai2=-4ki/m^ 

*’l2^^-52-^|)k^/m^ 

a21 = {?2-^l)Ci/Ij(j( 

b2i = (Jl^+^)c^/Ix^ 

a22=(^2'-?P'^l/^xx 

b22=«^^2)'<^/Ixx-^>’/I 

aji=2(b-a)c^/I^^ 


a32=2(b-a)ki/I^^ 



( 12 ) 

(13) 


(14) 


(15) 


12 


— c^1=2(a-b)c^/ni^ 

C. , . 

' .' C2i=2(a^+b^)c^/Izz i ' ' 

C32=2(aW)k,/I^^-w^h/I^^ 

Equations (15) are a set of second order linear differential equations 
to be solved simultaneously. A digital computer program titled 
DYNAMO (front DYNAmic MOdels) was used to solve equations (15). 

Appendix A is a listing of the program used. " ' - 

Various terrain conditions may be simulated on the mathematical 
model. A step input to the vehicle c.g. is simulated by assuming 
an initial condition on the vertical displacement. An initial condi- 
tion on pitch displacement simulates a step either to the front or 
rear wheel pair as would be encountered by the vehicle's front wheels 
" hitting a sudden difference in elevation. Similarly an initial 

condition on roll displacement simulates a step to either the 
right or left side wheel pair. A washboard terrain effect can be 
modeled by equating the left side of the first of equations (16) to 

’ . •• i 

a sinusoidal forcing function of desired frequency and amplitude. 

The c.g. response to various terrains may be simulated by equating 
the left side of the proper equation of equation (15) to the appro- 
priate function. 

: This mathematical model may easily be extended to: 

(1) include a spring and damper to be located in the connection 
of the frame to the front axle, or 

(2) to consider wheels of different elastic properties. 

Appendices B and C list the resulting equations of motion. This model 

also be extended to include additional degrees of freedom. 



Specifically, because of a lateral stiffness problem with the all 
metal toroidal wheel, the additional degrees of freedom of lateral 
translational and yaw displacements may be useful extensions. 



• CHAPTER 3 

. MEASUREMENT OF THE 0.4 RPI-MRV 
• PHYSICAL CHARACTERISTICS 

Once the equations of motion were developed, the physical 
characteristics of the 0.4 scale model MRV were determined in order 
to convert the equations of motion into meaningful results. Most of 
the measurements were determined in a straightforward manner. The 
sprung, unsprung, and total vehicle mass (earth-weight), and the 
dimensions locating the vehicle center of gravity are listed in Plate 2 
The remaining characteristics proved to be more difficult to measure, 
and are discussed below. 

Spring Constant . The sole source of flexibility considered in the 
mathematical model was the all-metal wheel. All four wheels were 
tested with the apparatus shown in Plate 3. The force-deflection 
curve of figure 5 shows the results. The spring proves to be quite 

I 

linear with a spring constant of k-=34.5 lb_/in. As discussed ' 

If- • , 

later, the actual dynamics of the MRV were measured with a capaci- 
tance-displacement measurement system, capable of measuring dis- 
placements of under one inch. To protect the capacitance probes, 
displacements were limited to a maximum amplitude of 1/4 inch. 

Under such small loads yielding small deflections, the spring con- 
stant of the wheel was found to be lower than the large-load spring 
constant. Figure 6 shows the wheel deflection characteristic 
under small loads, indicating a spring constant of k^=24.8 Ib^/in. 
Damping Coefficient . The damping coefficient of the wheel, associ- 



Summary of 0.4 Scale RPl-MRV Physical Characteristics 


Vehicle 

1 

Total earth weight 

w=51.0 Ib^ 

Total mass 

m-1.58 slugs 

Front wheel load 

7.5 Ib^/wheel 
r 

Rear wheel load 

18 Ib^/wheel 

Sprung earth weight 

w =49.0 lb. 
. r f 

Sprung mass 

m^=1.52 slugs 

Unsprung earth weight 

w =2.0 Ib- 
u f 

Unsprung mass 

m =0.06 slugs 
u ® 


Dimensions: (see Figures 1 and 9) 

Jt^^=21.0 in 
J^2=23.0 in 
a=22-l/4 in 
b=34-3/4 in 
h=10.0 in 

Moments of Inertia of rolling mass (ra^) (see Figure 2) 


v*^.=2.73 

slug-f t^ 

I =3.73 slug-ft^ 

XX 

XX 

: , ,=6.30 

slug-f t^ 

I =6.30 slug-ft^ 

y y 

yy 

: , ,=5.67 
z z 

slug-f t^ 

I =6.67 slug-ft^ 
zz 


Wheel 

Spring constant k^=414 Ib^/ft (large loads) 

k^=298 Ib^/ft (small loads, see discussion) 
Damping coefficient Cj^=4.0 Ib^/ft/sec 


PLATE 2 



/^orce (/6) 
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Wheel Spring and Damping Tester 


ted with the spring, was more difficult to determine. The appara- 
tus shown in plate 3 was used. Several tests were conducted in order 
to find a few tests during which the wheel was allowed to vibrate 
freely. The capacitance displacement measuring equipment was used 
to monitor the motion. A strip chart recorder was used to gain a 
time displacement plot (figure 7). In order to decrease the period 
of oscillation, and therefore increase accuracy, the tests were 
conducted with the wheel axle loaded to 8 pounds. The critical 
damping coefficient is easily determined from 


C^=2rajk7^ 

The logarithmic decrement (<T) of vibration decay is determined 
by measuring the amplitude of successive peaks of vibration (figure 7) 
and the following expression 

The logarithmic decrement can then be related to the damping ratio 


by 


=r- 


Since the d imping ratio is the ratio of actual to critical damping 

=C/C^) , the actual damping coefficient (labeled for the wheel) 

is then determined. The damping coefficient of the wheel was 

determined to be C =4.0 Ib^/ft/sec. 

if 

Moments of Inertia . The bifilar pendulum method was used to deter- 
mine the vehicle moments of inertia about its center of gravity. 
Since the sprung mass is much greater than the unsprung mass 

(m^) , the moments measured were taken to be the moments of the 


sprung mass of the vehicle. Plate 4 shows the arrangement used. 




FIGURE 7 


Damped Wheel Response 
Under 8 Pound Load 



Sc^=H.O 

C- H.O Ibf /•Pf'/seo 




I 

I 



of I for the 0.4 RPI-MRV 

I 

I 



The commonly used method Is basically a torsional pendulum. The 
vehicle must be suspended so that one of its axes is vertical. 
Suspended by two identical lines of length L attached to the vehicle 
at distances and from the c.g. (see Figure 8), the vehicle is 
displaced angularly in the horizontal plane about the vertical 
axis through its c.g. A natural torsional pendulum motion results. 
Once the period of oscillation is determined, the moment of inertia 
about the vertical axis through the c.g. is determined from 

where w is the total vehicle weight and t is the period of os- 
cillation. This measurement procedure was carried out for all 
three axes of the primed coordinate system of Figure 2, Plate 2 
tabulates the results. The moments of inertia about the unprimed 
coordinate system of Figure 2 are obtained from the primed axes 
by use of the parallel axis transfer expression 

where h is the distance separating the two axes (Figure 2). These.. 

values are also tabulated in Plate 2. 
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CHAPTER 4 

. DYNAMIC TESTING • ’ 

Approach . The purpose of the dynamic testing was two- fold: 

(1) Data was sought for a comparison with the predicted response 
by the mathematical model 

(2) The sources of undesirable motion were located in order . 
to indicate possible corrective alterations. 

The dynamic testing proceeded in two directions: 

(1) Measurement of transient responses to step input dis- 
placements 

(2) Measurement of frequency responses to forced constant 
amplitude sinusoidal displacements 

Responses were measured at various vehicle locations as a 
result of various inputs at different locations. It was then readily 
apparent that the testing procedure could be simplified. Measure- 

• * • i 

■ ! 

ment of the response of the c.g. provided all the necessary infer-* 
mation. The procedure is described below. ^ 

Equipment . All measurements were made with a capacitance displace- 
ment measurement system. Capacitance probes with a maximum range 
of one inch were used. Plate 5 shows a probe in place for testing. 
After the signal from the capacitance probe was amplified, it was 
connected to a strip chart recorder for an output plot. All dis- 
placements were calibrated with a dial indicator. 

An MB electromagnetic shaker unit was used to provide the sin- 

A 

usoidal forced displacements. The shaker unit was fitted with a 
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PLATE 5 


Close-up View of the 

Capacitance Displacement Measuring Probe 


wood extension arm to provide two platform pads on which the test 
vehicle was placed (Plate 6) . Initially the frequency response of 
the shaker unit with platform was measured in order to separate the 
vehicle response from the shaker response. As shown in Figure 9, 
the shaker unit frequency response is not particularly flat. How- 
ever, no sharp peaks exist. A comparison of the MRV frequency 
response (Figure 13) with the shaker unit frequency response shows 
that the shaker unit caused no gross errors or false peaks. 

Testing . 

Step Response . The transient response of the c.g. of the MRV 
to step inputs was measured, Plate 6 shows the general equipment 
setup for measurement of the vertical response. A close-up view of 
the capacitance probe and its location relative to a conducting 
surface on the MRV is shown in Plate 5. The roll and pitch responses 
were easily determined by the use of two capacitance probes. The 
two probes were placed on a line parallel to or coinciding with the 
principal axes of the vehicle. (Plate 7 shows the equipment setup 

and capacitance probe locations for pitch angle response measure- . 

( 

ments) , For small angles of rotation, the change in angle was 
simply the difference of the probe readings divided by the distance 
between the probes. The amplifier used to amplify the signals from 
the probes was equipped with a difference junction. Therefore the 
response was easily obtained by connecting the difference junction 
to the strip chart recorder for a time plot. 

All step inputs were simulated by the release of an initial dis- 
placement, Using nylon line, a weight was suspended from the desired 
location on the MRV, The nylon line was then burned to release the 
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PLATE 6 

Equipment Setup for Vertical Displacement 


Response Measurement 
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PLATE 7 


Equipmemt Setup for Pitch Angle 
Response Measurement 



weight* Figures 10,11, and 12 show step responses of vertical, pitch, 
and roll displacements respectively. 

Forced frequency response . The MRV may be traversing Martian 
terrain similar to a washboard surface. Therefore, it is important 
to determine the natural frequencies whiVh exist in the MRV, and the 
relative amplitude of vibration at these frequencies. As explained 
above, the electromagnetic shaker unit was used to force inputs of 
various frequencies. Since the shaker unit is of too small capacity 
to vibrate the 0.4 MRV as a whole, the vibration tests were conducted 
in two phases. Only one set of wheels, front or rear, was driven 
at one time. Figure 13 shows the frequency response to the sinusoidal 
inputs to the front and rear wheel pairs. Each curve shows two 
peaks at the frequencies of about 4.0 cps, and 6.5 cps. The location 
of the peaks in the response is easily predicted by evaluation of the 
expressions 

j -9- (17) 

If the total load on the rear wheels and the equivalent spring 
constant of the rear wheels (2k^) is substituted into equations (17), 
the result of f=4.34 cps is obtained. Likewise, if the total front 
wheel load and spring constant is used, f=6.74 cps. These predicted 
natural frequencies are reasonably close to the measured frequency 
response peaks of 4.0 and 6.5 cps. The association of the 4.0 cps 
peak with the rear wheel pair and the 6.5 peak with the front wheel 
pair is also a logical conclusion from Figure 13. In the case of the 
rear wheel pair input curve, the higher peak is associated with a 
frequency of 4.0 cps, which results from the rear wheel spring. The 
6.5 cps peak is similarly associated with the front wheel spring 
on the front wheel pair input curve. 




FIGURE 10 

Vertical Step Response — Actual 



Pitch Step Response— Actual 





FIGURE 12 

Roll Step Response — Actual 



0.4 Scale MRV Frequency j Response 
• • (vertical response of c.g.) 
Displaceraent vs. : Frequency : ; | j 











CHAPTER 5 

COMPARISON OF MATHEMATICAL MODEL WITH TEST RESULTS ■ 

AND CONCLUSIONS 

Figures 10, 11, and 12 show actual step input responses of 
the 0.4 RPI-MRV 1n the vertical, pitch angle, and roll angle modes, 
respectively. By setting the initial conditions of the mathematical model 
equal to the initial conditions imposed upon the physical model , a 
.mathematically predicted response was obtained. Figure 14 shows 
the predicted vertical response while Figure 15 shows the pitch and 
roll responses. 

A close comparison of the actual response with the predicted 
response indicates the following: 

(a) The period of oscillation of the actual model is about 
0.3 seconds in the vertical displacement and pitch response and 
0.33 seconds for the roll response. The corresponding values of 
the predicted response are 0.2 and 0.24 seconds respectively. The 
mathematical model predicts a somewhat faster response, but correctly 
predicts an approximate relative difference between the roll response 
and the vertical /pitch response. 

(b) The ratio of the first and second "peak" of the actual 

vertical response is mathematical model predicts 

less damping showing the same ratio to be 1=2.1 . However, 

the general envelop of the actual decaying response shows the oscil- 
lation to be reduced to 20% of the initial value after 0.6 seconds. 

The mathematical model predicts the oscillation envelop to be at 
about 15% of the initial peak value at the equivalent time. Also 
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.to be noted Is the Increase in amplitude of oscillation of the actual 
vertical response between the fourth and fifth peaks in Figure 10. 

This is due to the transfer of energy from the roll and pitch modes 
to the vertical mode. The mathematical model correctly predicts the 
increase in amplitude between the fourth and fifth peaks. A compari- 
son of the actual versus predicted roll and pitch angle responses 
yields, similar results. The general envelop of oscillation is similar 
in each case. However, the increase in amplitude of the actual re- 
sponse, because of the energy transfer, is not reflected in the pre- 
dicted response. 

Sources of error in the mathematical model are most likely due 
to the follov-zing mathematical deficiencies: 

(1) The mathematical model is based on a point mass located 
at the c.g. of the sprung vehicle. The actual model is quite dif- 
ferent with no one general area of mass concentration. 

(2) Only three degrees of freedom have been considered in the 
mathenatical model. In reality, the degrees of freedom not considered 
(lateral and longitudinal horizontal translation and yaw angle dis.- 
placement) do experience considerable excitation. Since these modes 
are coupled with the three degrees of freedom considered, there 

will be considerable energy transfer among all modes. Therefore, 
neglecting three modes induces considerable error into the mathema- 
tical, model. 

(3) Many of the members of the 0.4 scale model are rather 
loosely joined. This "slop" is obviously not included in the mathe- 
matical model and is therefore a source of error. 

An overall evaluation of the mathematical model shows it to 



be of some value in predicting gross changes in the vehicle response 
as a result of design parameter changes. For instance, should it 
be desired to apply a wheel of different spring and damping coef- 
ficients, the mathematical model may be used to predict general 
responses from which comparisons may be made. Likewise, changes in 
wheel base dimensions or payload location may easily be translated 
into the mathematical model. 

In the optimization of the final vehicle design, a more accurate 
mathematical model will most likely be desired. The model pro- 
posed herein will provide a good basis for extension. Additional 
degrees of freedom could be included to provide a more sophis- 
ticated and more accurate model. For the present time it is felt 
that this mathematical model will be adequate to help in determining 
the main design parameters of the 0.4 scale physical working-model, 
the construction of which is the next step in RPI's vehicle develop- 
ment program. 
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APPENDIX A 


SAMPLE DYNAMO PROGRAM FOR EQUATIONS (15) 


11/05/73 

I Yl.K-Yl.J+nT*YlDOT.JK 
R YlU0T.ia = Y2.K 
L Y2. K=Y2. J+DT*Y2D0T, JK 

R Y2D0T. KL=/M1*Y2. K+A12*Y1. K+B11*R2. K+B12*R1. K+C11*P2. K+C12*P1. K 
L Rl.K = IU.J + t)T*fUDOT.JK . 

^ JR RIDOT. KL = R2. K ' 

R2.K = R2. J + DT*P^2D0T.JK 

R R2DOT.KL = A21*Y2.i;+A22*Yl,K+B21*R2.K+B22*Rl.K 
L P1.K=P1. J+DT*P1DUT. JK 
W P1UUT,KL = P2. PC 
L P2. K=P2. J + nT*P2D0T. JPC 

.R P2D0T. 1CL = A31*Y2. K + A32*Y1, K+C31*P2. K + C32*P1JC 
fj All--4.0*Cl/IiR 
N A12--4. fJ*ICl/flR 
n Ull = ( 1.1' L2)*C1/MR 
n B12 = (Ll-'L2)-MCl/r!R 
M Cll = 2. 0*(A"B)*Cl/riR 
N C12 = 2. P*(A-B)*ia/n[i 
n > II a21 = (L1-L2)*C1/IXX 
^ ] ;j A22 = (L1-'L2)*K1/IXX 

U n2i = -(Ll^-Ll + L2*L2)*Cl/IXX 
U B2 2 = ~( ( Ll^ai + L2*L2)*fCl/ I XX) + (UI’*H/ I XX ) 

U A31-(A-R)*2.0*C1/IXZ 

U /.52 = (A-C)*2.0*lCl/!/-Z ( 

n C31 = -(A*A + [UB)*2. 0*01/ IZZ ■■ f 

M C32-~((A*A+B*B)*2.0*ia/IZZ)+a/R*H/IZZ) ' ■ 

' C 01-02. 0 ' ■ ' 

C Kl-290.0 
0 MR- 1.52 
N A-22. 25/12.0 
IJ B-34. 75/12.0 
N Ll-21. 0/12.0 
^ I U L2 -2 3. 0/12,0 
. C lXX-3.73 
C IZZ-0.C7 
C Uli-40.0 
, N M-10. 0/12. 0 
' U Yl-O.OUIG 
n Y2-0.0 

p. u ui- 0 . oior> 

^ I fJ IJ2-0.0 - 
H PI -0.0293 
^ U P2-0. 0 
PRINT Y1,R1,P1 
£ < PLOT Yl-Y 

PLOT R1-R,P1=P 

SPEC UT-0. 001/ LENGTH-01. O/PRTPFR-0.0 2/ PLTPER-0. 02 
RUN 
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• j ■ . 

' ■ I 

■ Brief Explanation of DYNAMO Program 
Statements . , ' 

A The set of three second order differential equations 

(equations (15)) expressed as six first order dif- 
ference equations in the form required for DYNAMO. 

B . The constants of equation (16). 

C The physical characteristics of the 0.4 scale 

RPI-MRV. 

D The initial conditions. 

C The output statements. 









•APPENDIX B 


ALTERNATE MATHEMATICAL MODEL— 1 

To increase the damping and to improve the overall vehicle 
response, it has been proposed to join the MRV frame to the front 
axle through a spring and damper. Figure 15 shows the schematic 
diagram for such a configuration. The following is a derivation 
of the equations of motion for this model. 

Inertia Forces . The inertia forces are the same as derived in 
Chapter 2. 

External Forces . From Figure 17, the equivalent rear and front 
mechanical impedances respectively are: 

.2’^= 3(fc, 4C|s) 

_ j(ki-»C|s)(l<a»Cjs) 

■- ' P Uc,+Ca)s 

Following the same procedure of virtual displacements as developed 
previously: 

(a) roll with vertical displacement (Figure 3) 

~3A,A^y (19) 

(20) 

where A ~ 

4 = 

As= a(k, + c,s) 

(ki, + Cjs) 

m 





(a) Rear suspension 



2,= aOc.^c.s) 

<^2 = k, - CjS 
z.f = (^k|-» 3C[a)(l<j'*' yCj) 

(b) Front suspension 

FIGURE U 

Alternate Model of RPI-MRV 
Front and Rear Suspension 
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(b) pitch with vertical displaceinent (Figure 4) 

( 21 ) 

y (b^r^-aZ^)*e^(-b^2^-a’z^ + hWr) (22) 

Combining equations (5), (19), (20), (21), and (22) and simplifying, 
the equations of motion take the following form: 


§»•••* ** * ** * * 

^ • • •• • 

Cj, 4 1 c3^e^*-CssQ^* Oj, y + y f) 3 i ( Sx) •* ®x) = 0 


where: 

a„- -[z^.][(3k,vk,) + -i;,(Mc,c^^pc,=)j 

Q(a = ■■ [cJTi^,] [ ( l^)(lc,Cj4 c, kj tc, 0] 

*^13 ”~Lcj*'sC,]{( 

bn “ * [cj'tsc,] [ ( "mr) 

^la ~ ’*'Lc^->s?,l[k('"j‘*^‘-i)+C|(3k|-i (r^^( w,.)(-?a~'?/) j 

t>i 3 " [,cj 1 ac,l ( ^ J '-Pi) ("S !</ + fei ks} ' 

~ ■* Lcj\sc,l [pbC|^ + “b)c,Cj] 

C.a - + [-£^,] (l^j [Mc,b k, (<^- b)Cc^k,+ c, Q] 

^13 ~ lcj,-tx7] ( m)) L^i ki (a - b) + Si> k,^j 

Qa. = L^J ( (-Pa '-PJ (^- O 

Qaa " + ( i;^)(^>a--P.) ( b,Cj 4 Jr^c,) 

<^a3 ■ ■‘Lc^^acji (_Xjj,^)('Pa'-P<)(kka) 

t>ai =^[^io^;][(i?k,+ U+ ^(■^i^-Pj'Xc.Ca + ^C,^)] 

t>ja ■ "^Lcj^ac,] [(-Pi'^^ia )[Ci(?k,-^kj)+k, ( c4»3c,)] - (c/f^, i ej)| 
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0^1^ ■** I C;i*f 5c,l ( 

^(b-o)(C|k;j + C^k,)J 

^33 ^ [ c^+acj ( Xx^s. ) +^kika (b'o)] 

Cjj ~ Lc;j-»;?c, j |(^k,^kj)+ [b(o 

^3:t ~ *^ E C;j 4 5C|] ( ^ oka+C^ki) *" iVf^h f3c,-tC^^ 

<^J3 ^ L [Ha^t>i<'i^'+ 2>Q^fc>t|kj‘^4tfef +^b% fea 

"- W|- b ( ^ t'l ■*’ )] 

Klji = - [ c^^,] ( ) {^*=1+ ka) 

^31 ' [ (“ (<3k,+ fcii) 

n^3 = ~ ^Txy/Xzz. 

The digital computer program DYNAMO is an excellent method to use 
to solve equations (23). 

Initially one' would believe that the simpler mathematical model 
developed for the schematic of Figure 1 could be obtained from 
equations (23) merely by letting' the values of and approach 
infinity. A look at the coefficient terms shows that no meaningful 
expressions result if this limit is taken. An alternate way to 
negate the effects of and ^2 would be to equate them to high 
values in the solution of equations (23) with the DYNAMO program. 



Although this solution would yield desired results, it would be 
impractical. Large values of k and C associated with relatively 
small values of m in mechanical systems yield high time constants 
(time constant T «. 1/4 17m) « Any method of solution of differential 
equations depends upon a sufficiently small step (time interval 
in this case) to yield reasonably accurate results. A smaller step 
size necessitates a larger number of iterations, in the calculations 
to span a specified time interval. This results in considerable 
more computer time. Therefore, equations (23) should not be used in 
order to model the simplified schematic of Figure 1. 



• • APPENDIX C 

. ■ . I ■ 

■ ■ . . ALTERNATE MATHEMATICAL MODEL— 2 

Because of limitations at the time of construction of the 0.4 
scale MRV, all four toroidal wheels have identical physical properties. 
The mathematical model developed in the main text of this report is 
based upon this existing model . ‘ — 

Since the vehicle configuration is designed so that the rear 
wheel pair supports 70-80% of the total vehicle v^eight, the stiffness 
and damping coefficient of the rear wheel should be proportionally 
greater than the front wheel properties. This appendix presents a 
mathematical model where the physical characteristics of the wheels 
are not identical. The subscript ‘f* denotes front while the sub- 
script 'r' denotes rear. 

The inertia forces of equation (5) remain the same. The ex- 
ternal forces due to displacement in the vertical and roTl mode are ^ 
rewritten from equations (6) and (7) as follows: . I 

fy = "(kr»Crs)(y--PA)-(l<c-*Ci-s)(y+J?a©x)' 2(k| * C^s) 

+ K(h-y)&^ 

Using the small displacement approximations, JJ •j» h0^, 

and h»y, equations (24) and (25) become 

( fe r Cr S) [y (jp, -J),) - ©X ^ 


(27) 
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Similarly the resultant external forces due to a displacement In the 
vertical and pitch mode follow from equations (10) and (11): 

Fy= -a(krrsC^)(y-QeJ-5(k^tC^s)(y + tej;) (28) 

^rs)(a~h&^){Y-a&2)~3Q:ftCfs)(b+h&£)(ytb6^) 

+ VVr(h-y)^ (29) 

Again with the small angle approximations these equations reduce to 

Fy= “<^[_(kr‘»-Cr5)(y-Q0£)4^(k^*»C^s)(y + (30) 

l^z.- 3ft(MCrS)(y-aei)-5b(k^+Crs)(y +t>0^) +vv^h0^ (31) 


Combining equations (5), 

(26), (27), (30), and (31), and 

rearranging 

the final form results: 





y 4 o„y +o,^y 


C„©Z + C,J0^ 

= o 


©x'*’ kii^K + 

+ Qj, y •» 

Q^jy=o 




*Oj,y =0 


where 







K ■■ 




~ 1iv 


"U''*‘’/jkr/rxx-W, 

^/Xxx 

hi = 








■ tf/Xxx 



^ (q "I^) ^r/ fvjf- 



C IP = 

3(Q-t) kr/mr 


(Q’C^ + fc*C^) 






w,h) 









a 

" Cbfef.-aSr^) 
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k^=sprlng rate of rear wheel 
k^=spring rate of front wheel 

C^“damping coefficient of rear wheel . . 

C^=damping coefficient of front wheel 

a,b,fj ,Jf 2 jh, refer to the dimensions of Figure 1 

refer, to vehicle physical characteristics of 

Plate 2 
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